DNA Replication via Entanglement Swapping 



Onur Pusuluk 1 and Cemsinan Deliduman 2 ' 3 

1 Istanbul Technical University, Faculty of Science and Letters, 
Physics Department, Maslak 344^9, Istanbul, Turkey 
2 Mimar Sinan Fine Arts University, Department of Physics, Be§ikta§ 34349, Istanbul, Turkey and 
3 Feza Giirsey Institute, Qengelkdy 34684, Istanbul, Turkey 

Quantum effects are mainly used for the determination of molecular shapes in molecular biol- 
ogy, but quantum information theory may be a more useful tool to understand the physics of life. 
Organic molecules and quantum circuits/protocols can be considered as hardware and software of 
living systems that are co-optimized during evolution. We try to model DNA replication in this 
sense as a multi-body entanglement swapping with a reliable qubit representation of the nucleotides. 
In our model molecular recognition of a nucleotide triggers an intrabase entanglement corresponding 
to a superposition state of different tautomer forms. Then, base pairing occurs by swapping intra- 
base entanglements with interbase entanglements. We examine possible realizations of quantum 
circuits to be used to obtain intrabase entanglement and swapping protocols to be employed to ob- 
tain interbase entanglement. Finally, we discuss possible ways for computational and experimental 
verification of the model. 

PACS numbers: 03.67 Ac, 82.39.Jn, 87.14.gf, 87.14.gk 



According to the central dogma of molecular biology, 
genetic information stored in DNA is duplicated by repli- 
cation and is used by successive transcription and trans- 
lation. During replication, enzyme DNA polymerase 
(DNApol) recognizes the nucleotide bases N = {A, T, G, 
C} of template DNA strand and finds their complemen- 
taries {A=T, T=A, G=C, C=G} from the surrounding 
environment for base pairings. Recognition interaction 
between this single-stranded DNA (ssDNA) and poly- 
merase enzyme is one of the several unknown aspects of 
replication. Also, mechanism used for finding the correct 
nucleotide from the surrounding environment is a mys- 
tery. Since a lot of amino acids exist in the active site 
of DNApo/ pQ, both experiments and quantum chemi- 
cal calculations are insufficient to clarify these mysteries. 
However, there are some quantum information processing 
models proposed for replication mechanism. 

In 2001, Patel [2] formulated nucleotide selection from 
surrounding environment as an unsorted database search 
and examined the possibility of the use of Grover's al- 
gorithm [3 in the evolutionary context. Although he 
achieved to model base pairing as oracle in the algo- 
rithm, initiation of the algorithm requires the superpo- 
sition of four nucleotides which is not quite possible. In 
the wave analogue of the algorithm [4J replication should 
begin with the interaction of DNApo/ and all of four nu- 
cleotides, but it is known that DNApol first binds to 
DNA template in this process [T]. Recently, Cooper [5] 
showed that molecular genetic transcription data of bac- 
teriophage T4 is compatible with a quantum treatment in 
which enzyme makes a measurement on coherent protons 
that causes an entanglement between enzyme and pro- 
tons. According to this model, enzyme-nucleotide inter- 
action takes place on the protons and electron lone pairs 
that contribute to the base pairing. However, Watson- 
Crick (WC) edge of the nucleotides (Figure [I]) can not 



be represented as four orthogonal states by such a treat- 
ment. In this paper, we try to model DNA replication 
as a multiparticle entanglement swapping [6 with a re- 
liable qubit representation of nucleotides. In the model, 
molecular recognition of a nucleotide takes place on the 
protons and electron lone pairs that are present on the 
Hoogsteen (H) edge. 

Meanwhile, covalent structure of the nucleotide bases 
is not static. Random, reversible and infrequent proton- 
coupled electron derealizations convert bases into rare 
tautomer forms (N*, G**, 0**). In 1963, Lowdin [7] 
claimed that tautomeric shifts on both paired DNA bases 
by a double proton tunneling through hydrogen bonds 
can be responsible for mutations. After some corre- 
lated ab inito calculations [5J |9] that support Lowdin's 
claim, Villani [lOj [IT] calculated more reliable poten- 
tial energy surfaces by density functional theory (DFT) 
method and found possible AT^A*T*, G-C^G^C* 
and GC— ^G*C* transitions by concerted or two step 
double proton-coupled electron transfers. Since cal- 
culated transition probabilities are pure quantum me- 
chanical ones, tautomerization in double-stranded DNA 
should have a quantum mechanical nature. However, cel- 
lular environment may cause a decoherence effect on nu- 
cleotide base in both free nucleotide and ssDNA cases. 
Therefore, we assume that enzyme-nucleotide base in- 
teraction can suppress the decoherence effect and bring 
the state of nucleotide base to a superposition of all tau- 
tomer forms. Indeed, superposition of all tautomer forms 
is nothing else then intrabase entanglement of the atoms 
on WC edge. 

Base pairing occurs via hydrogen bonding, but nature 
of the interbase hydrogen bonds is not well understood 
because of the insensitivity of experiments. A nonlo- 
cal DFT method [12] showed that covalent contribution 
to hydrogen bonds is 38% in AT pairing and is 35% 
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in G-C pairing. This conclusion was then supported 
by subsequent DFT studies [13] and similar conclusions 
were reached by semiemprical methods with geometrical 
and atoms-in- molecules topological parameters, natural 
bond orbital analysis, and spectroscopic measurements 
[T4] . Covalency of hydrogen bonds means that electron 
lone pair of the acceptor orbital is quantum mechanically 
shared between its own orbital and unoccupied antibond- 
ing orbital linking donor and hydrogen atoms. Thus, we 
can interpret the state of a hydrogen bonded atom pair 
as an entangled state. Hence, in our model, bases are 
paired by swapping intrabase entanglements with inter- 
base entanglements. 
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FIG. 1: Parts of the nucleotides: atoms can be grouped 
according to region they will be found in DNA. Hoogsteen 
(<H>major groove), Watson-Crick (^pairing plane), and Sugar 
(<H>minor grove) edges are indicated respectively by dashed, 
filled, and plain boxes. Arrows inside the boxes show the 
order of the qubits used in qubit representation. 

We assume that base recognition by DNApol should 
occur over the H edges of the nucleotides via a quan- 
tum measurement. In consensus, hydrogen bond donor 
and acceptor atoms of bases are only O and N atoms. 
However, there are a small number of computational ob- 
servations in which C atoms of nucleotide bases have the 
ability to make blue-shifting hydrogen bonds [9]. In this 
respect, when electronic configurations of the individual 
O, N, and C atoms in H, WC, and S edges are consid- 
ered, it is found that each atom has two different energy 
states: a relatively lower energy state for acceptor situa- 
tion and a relatively higher energy state for donor situa- 
tion (Figure [2]). Surprisingly, only H of the nucleotides 
can be represented as four orthogonal states if lower en- 
ergy states are regarded as |0) qubit and higher energy 
states are regarded as |1) qubit: 



\A)h = \01), 
\G)h = |00> , 



\T) H = m, 

|C>» = |11>. 



(1) 



The first qubit carries information about purine- 
pyrimidine distinction, whereas the second one carries 
information about imino-enol distinction. In this sense, 
DNApol should pair bases whose qubit representations 
are complementary to each other. Not only correct base 
pairings, but also mispairings like AC* and G*T pair- 
ings can be accounted for by this assumption. 
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FIG. 2: Electronic configurations and qubit representations 
of the O, N, and C atoms: configuration indicated by * is 
not present in any tautomer form. However, it is possible to 
observe it in blue-shifting hydrogen bonds of DNA. 



However, since pairing occurs between WC edges of 
nucleotides, information processing starting with recog- 
nition over H edges should continue over the WC edges. 
Equality of the second qubit of H edge and the first qubit 
of WC edge makes such a transition in interaction re- 
gion reasonable. Qubit representation of WC edges of 
nucleotides before replication are as follows: 



|A) W c,/ = |101), 
\G)wcj = |011) , 



\T)wc,i 
\C)wc,i 



1010) , 
1 100). 



(2) 



A double proton transfer between the DNApo/ and the 
first atom of WC edge which occurs immediately after 
the recognition, can trigger a tautomeric transition (Fig- 
ure [3]). If such a transfer has a quantum nature, recogni- 
tion interaction can trigger an unitary transition to the 
superposition of all tautomer forms. A candidate for such 
a transformation U is shown in the Figure [4} Since |0) 
and |1) states of an atom respectively correspond to the 
absence and presence of a proton bonded to that atom, 
NOT gate can be regarded as a quantum mechanical pro- 
ton transfer. Due to the same reason, SP gates can be 
considered as formation of a quantum mechanical hydro- 
gen bond between the enzyme and the particular atom. 
Both the proton transfer and hydrogen bonding are the 
usual tasks done by enzymes and there are some evi- 
dences for the unignorable role of quantum effects and 
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FIG. 3: Tautomeric transition by proton transfer between 
enzyme and nucleotide 



dynamics on the enzymatic reactions [15 j. This means 
that intr abase entanglement by transformation U is a 
possible action for the enzyme DNApo/. 
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FIG. 4: A possible quantum circuit for U which transforms 
\N)wc,i states to the \N) wc.q states: superposition matrix 
SP(0) of controlled — Superposition gates equals to the mul- 
tiplication of rotation matrix R(9) and Pauli-Z matrix 

To provide an equilibrium between maximal entan- 
glement and robustness, we take the angles and <p 
in the quantum circuit for U as arccos(v / 2/v / 3) and 
arccos(l/v // 2), respectively. Then, \N)wc,Q states are 
obtained as: 



\A)wc,q = (|01> - 

\T)wc,q = (|01)^ 

\G)wc,q = (|011) 

\C)wc,q = (|100) 



|10»|1)/V2, 

|10»|0)/V2, 
f |101) + |110»/V3, 
- |010) + |001))/V3. 



(3) 



To consider each base pair as an intact system, tensor 
products of these states should be taken. However, we 



reorder qubits of these product states in such a way that 
hydrogen bonded atom pairs come next to each other in 
order to clarify base pairing. Then, we get 
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\A • T) wc ,q = 



ii) 

|10)|01) 



\G- CVc.q = -(|01) 1 10) 

-|00)|11) 
+|00)|10) 



|10)+|01)|10) 

|10)-|11)|00) 

|10)+|11)|00) 

|10)-|10)|01) 
|11)+|10)|00) 



|10> (4) 

|io», 

|10) + |11>|10>|00) 

|10> — 110>|11>|00> 
|11) + |10)|10>|01» 



If intrabase hydrogen bonds have a quantum nature 
as discussed, each hydrogen bonded atom pair of two 
paired nucleotides (N r N 2 )wc,o can be considered to 
be in an entangled state. Then, in order to turn in- 
trabase entanglements into interbase entanglements, U 
should be followed by an irreversible transformation, 
S : \Ni • N2)wc,Q \Ni • N2)wc,o, which is an entan- 
glement swapping (Figure [5]). We observe that, in the 
case of G-C pair, before S there are two three-qubit (in- 
trabase) entanglements and after S there are three two- 
qubit (interbase) entanglements. Similarly, in the case of 
AT pair, before S there are two two-qubit entanglements 
and after S there are two two-qubit entanglements. 




FIG. 5: Entanglement swapping model of replication 

Swapping intrabase entanglements to interbase entan- 
glements can be achieved by a five-step protocol S as 
follows: 

1. Third and fifth qubits of the reordered base pair 
states (second and third qubits of the nucleotide base in 
template DNA) are subjected to a transformation V as 
shown in Figure [6| 

2. A Bell measurement is performed on the third and 
fourth qubits of the reordered states. 

3. If the result of the measurement is one of the two 
Bell states |/? o> and ((|00) ± |ll))/\/2), fourth and 
fifth qubits of the reordered states are subjected to the 
Pauli-X transformation. 

4. A Bell measurement is performed on the first and 
second qubits of the reordered states. 
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5. If the result of the measurement is one of the two 
Bell states |/3qo) and |/?io), second and fifth qubits of the 
reordered states are subjected to the Pauli-X transfor- 
mation. 
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FIG. 6: Transformation V in the protocol which swaps intra- 
base entanglements to interbase entanglements: this transfor- 
mation entangles the qubits on the condition of their equality. 
Hadamard matrix H equals to SP(tt/4). 

If possible |7Vi • N 2 ) 

wco states are written in terms 
of the Bell states |An) and \j3 u ) ((|00) ± \ll))/V2), then 
state ensembles of base pairs after S are found to be 



|4-TVc,o = {0.43,|A)i>|A)i>|10) 
0.07,|/?ii>IA,i)|10) 

0.07,|/3 i>|/3ii)|10) 



\G ■ C)wc,o 



{Pj m ,\(3 jk m 



jm 



(5) 

O.43,|0ii)|/3n>|lO)}, 
■01) ' 
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TABLE I: Probabilities and probability amplitudes in state 
ensemble of G-C pair 



I 




^00 


pi 




&2,10 


pi 


1 


+0.51 


-0.86 


0.11 


-0.51 


+0.86 


0.11 


2 


+0.38 


+0.92 


0.09 


-0.38 


-0.92 


0.09 


3 


+0.96 


-0.28 


0.03 


-0.96 


+0.28 


0.03 


4 


-0.92 


+0.38 


0.02 


+0.92 


-0.38 


0.02 


I 


abi 


b l oi 


pL 






Hi 


1 


+0.51 


+0.86 


0.11 


+0.51 


+0.86 


0.11 


2 


+0.38 


-0.92 


0.09 


+0.38 


-0.92 


0.09 


3 


-0.96 


-0.28 


0.03 


-0.96 


-0.28 


0.03 


4 


+0.92 


+0.38 


0.02 


+0.92 


+0.38 


0.02 



We note that, besides U, S also consists of only pro- 
ton transfer and hydrogen bonding processes, excluding 
Bell measurements. Bell measurements can be thought 
as formation of a quantum mechanical hydrogen bond be- 
tween measured atom pair if the outcome state is |/?oi) or 
|/?n). However, when the state of an atom pair collapses 
to |/?oo) or |/?io), atom pair and DNApo/ can not sepa- 
rate from each other since total proton number of base 
pair does not remain constant after the measurement. 
Thus, Bell measurements should be treated as formation 
of quantum mechanical hydrogen bonds between the en- 
zyme and measured atom pair if the outcome state is 
|/?oo) or |/?io). In such circumstances, conditional Pauli- 
X transformations can fix the total number of protons 
on base pair and make atom pair - enzyme complex sep- 
arable. 



Neither U nor S is unique for the given model. How- 
ever, this is not a disadvantage since there are several 
DNApo/ species and families with different replication fi- 
delities. This diversity in replication fidelity of DNApol 
can be accomplished by different U and S pairs. 

Since all of the states in both U and S can be ex- 
pressed as proton transfer and hydrogen bonding, our 
model could be tested by repeating the scenario with 
quantum chemical computations. Experimental verifi- 
cation is also possible. Evolution of the \Ni • N2)wc,o 
states in the presence of double well potentials can be 
prevented by sufficiently decreasing the time period be- 
tween two successive pairings. Then, probability of point 
mutations due to the formation of rare tautomer forms 
may increase according to Equation [5] and Table [TJ 

Entanglement swapping may be a basic tool used by 
enzymes and proteins in the cellular environment. If so, 
similar models may be developed for amino acid - tRNA, 
aminoacyl-tRNA - mRNA, and amino acid - amino acid 
interactions in the protein synthesis. If successful mod- 
els for these interactions can be developed, then we can 
achieve a deeper understanding of the role of the quan- 
tum effects and dynamics on the cellular information pro- 
cessing. 
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